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Abstract 

A detection at the Large Hadron Collider of a light Higgs pseudoscalar would, if in¬ 
terpreted in a supersymmetric framework, be a smoking gun signature of non-minimal 
supersymmetry. In this work in the framework of the Next-to-Minimal Supersymmet¬ 
ric Standard Model we focus on vector boson fusion and Higgs-strahlung production of 
heavier scalars that subsequently decay into pairs of light pseudoscalars. We demon¬ 
strate that although these channels have in general very limited reach, they are viable 
for the detection of light pseudoscalars in some parts of parameter space and can serve 
as an important complementary probe to the dominant gluon-fusion production mode. 

We also demonstrate that in a Higgsfactory these channels may reach sensitivities 
comparable to or even exceeding the gluon fusion channels at the LHC, thus possibly 
rendering this our best option to discover a light pseudoscalar. It is also worth men¬ 
tioning that for the singlet dominated scalar this may be the only way to measure its 
couplings to gauge bosons. Especially promising are channels where the initial scalar is 
radiated off a W as these events have relatively high rates and provide substantial back¬ 
ground suppression due to leptons from the W. We identify three benchmark points 
that well represent the above scenarios. Assuming that the masses of the scalars and 
pseudoscalars are already measured in the gluon-fusion channel, the event kinematics 
can be further constrained, hence significantly improving detection prospects. This is 
especially important in the Higgs-strahlung channels with rather heavy scalars, and 
results in possible detection at 200/fb for the most favoured parts of the parameter 
space. 


1 Introduction 

The presence of an extra singlet superfield in the Next-to-Minimal Supersymmetric Stan¬ 
dard Model (NMSSM) (see, e.g., |T] for a review) as compared to the MSSM, has a significant 
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impact on the ensuing phenomenology of the Higgs sector at the Large Hadron Collider 
(LHC). In particular, the NMSSM Higgs sector is enlarged by two neutral mass eigenstates, 
one scalar and one pseudoscalar, in addition to the three MSSM-like onesQ The singlet 
nature of the scalar component of the additional superheld allows Hi and Ai to be very 
light when singlet dominated, even down to a few GeV, without entering into conflict with 
current theoretical and experimental constraints. This is so because their couplings to the 
fermions and gauge bosons of the SM are typically much smaller than those of the doublet- 
dominated Higgs bosons (^^ 2,3 and A 2 ), which are assumed heavier. As a consequence, 
the observation of any of these potentially light states, in addition to the SM-like Higgs 
boson already discovered at the LHC Eia, would constitute a hallmark signature of a non- 
minimal nature of supersymmetry (SUSY). Careful examination of their mass and coupling 
values in relation to the mass and coupling measurements of the 125 GeV SM-like Higgs 
boson (and possibly other discovered Higgs states) in a variety of production and decay 
channels will eventually enable one to prohle their possible NMSSM nature. 

The lightest pseudoscalar Ai, in particular, can be the most crucial probe of the NMSSM 
as it can be very light, so that in principle it is accessible in meson decays, where it has been 
searched for initially laiaiiiiziEiisi- The Ai state with mass > 5 GeV has also been probed 
in the possible decay of a heavy (SM-like or not) scalar Higgs boson at LEP2 [lOl ITT] and 
Tevatron ca, where no signihcant excess was observed. Regarding the LHC, the situation 
is as follows. CMS searched for a light pseudoscalar produced either singly |13j or in pairs 
from the decays of a non-SM-like Higgs boson |TT] and decaying into the channel, 

while the AiAi —>■ 4 t signature (via a SM-like Higgs decay) has been tackled in [I5]j3 For 
completeness, let us also mention an ATLAS search for scalar particles decaying via narrow 
resonances into two photons in the mass range above 65 GeV |18) . though this does not 
constrain light pseudoscalars very much. 

In addition, there are plenty of phenomenological analyses aiming at assessing the scope 
of Ai discovery within the NMSSM at the LHG. Prior to the SM-like Higgs boson discovery, 
quite some effort was put into extending the so-called ‘no-lose theorem’ of the MSSM — 
stating that at least one Higgs boson of the MSSM would have been discovered at the 
LHG via the usual SM-like production and decay channels throughout the entire MSSM 
parameter space |19| (see |2nj for a recent review) to the case of the NMSSM |21112^ l23l 1241 
[25l [26l ET] [28l EHl 13011311 [ 321 13311311 [351 [36] . In the light of the recent Higgs boson discovery 
though, the above theorem is necessarily verihed and, if one wants to prove the NMSSM to 
be a viable alternative to the MSSM, one ought to probe it away from the MSSM limit. 

^Hereafter, our book-keeping of the physical Higgs states of the NMSSM will be as follows: the CP-even 
states will be denoted by Hi (with i = 1,2,3 and such that ruHi < < mnA while the CP-odd ones 

by Ai (with i = 1,2 and such that rriAi < mA^)- We also use Hsm for the scalar playing the role of the 
discovered Higgs boson and Hs for the singlet dominated scalar. 

■^The sensitivity of the di-photon sample to a singly produced Higgs boson (of mass 150 GeV and above) 
decaying to 77 has also been investigated by CMS 1161 alongside that of the 46 sample to pair production 
of 4 Lv Higgs bosons HZl. « ^ ^ 
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Following this line of reasoning, if one abandons the limiting case of SM-like decay 
channels of Higgs states, the NMSSM offers interesting signatures which are precluded 
in the MSSM after the latest experimental constraints, in the form of a variety of Higgs 
—Two-Higgs and Higgs —)• Gauge-Higgs decays. A large volume of phenomenological 
literature exists on these topologies, claiming that, for certain NMSSM parameter choices, 
these would be accessible at the LHC, eventually enabling one to disentangle the NMSSM 
from the MSSM hypothesis, thereby establishing a so-to-say ‘more-to-gain’ theorem |31| . 
The importance of such decays in the context of the NMSSM has been emphasised in 
Refs. |371 [38l |39] from the point of view of fine-tuning as well as a distinctive NMSSM 
signature at colliders. In particular, the Hi ^2 AiAi mode has received much attention. 
This decay can in fact be dominant in large regions of the NMSSM parameter space. It was 
realised that vector boson fusion (VBF) could be a viable production channel to detect the 
above modes at the LHC, in which the CP-odd Higgs pair decays into jjT~^T~ - [MlEailQ] 
(where j represents a jet). Some scope could also be afforded by a 4r signature in both 
VBF and Higgs-strahlung (HS) off-gauge bosons |28| . The gluon-fusion (GF) channel too 
could be a means of accessing Hi —?• AiAi decays, as long as the light CP-odd Higgs 
states both decay into four muons |30) or two muons and two taus m Finally, the 
scope of NMSSM neutral Higgs boson production in association with bh pairs was assessed 
in |22l [33l [23l 1421 I43l I44j . including the case of H 2 —)• ZAi decays with Z —)• jj and 

Ai —)■ T^T~ . 

All the above mentioned analyses were, however, performed prior to the discovery of 
the Higgs boson at the LHC mm- In the aftermath of the discovery, detection prospects 
of all NMSSM Higgs bosons, including also via their decays into other Higgs states, were 
recently investigated in |45j . though limited to the case of inclusive rates. Further, in |46j . 
the Ai —?• 77 decay channel was studied in the regime of a light Ai. In m it was then 
noted that in the NMSSM the Ai could in fact be degenerate in mass with the SM-Iike 
Higgs boson Hsm- It could thus cause an enhancement in the Higgs boson signal rates near 
125 GeV in the 77 , bb and channels simultaneously, provided that it is produced 

in association with a bb pair. The bbAi channel was also studied in detail in [48], with a 
more optimistic conclusion as compared to |4^. In |5^ the H^m AiAi —?■ Ai (with i 
denoting and fj,^) process at the LHC was studied in detail, while the bb^fi final state 
was deemed the most promising in |51j . The production of Ai via neutralino decays has 
also been recently revisited in |52[l53l lM]. Finally, NMSSM benchmark proposals capturing 
much of this phenomenology exist in the literature |35[ [55] . 

In this work we continue to pursue a recently started endeavour |49], with an intention 
to systematically analyse all the production and decay processes that could potentially 
lead to the detection of a light NMSSM Ai at the LHC with ^/s = 14TeV. In |49| (see 
also [MIET]), we considered the case of the light pseudoscalar Ai produced from a heavy 
Higgs boson coming from GF which then decays into either Ai pairs or ZAi (with the Z 
in turn decaying into electrons/muons). We found that the Ai can be accessible through a 
variety of signatures proceeding via Ai —>■ r+T“ and/or bb, the former assuming hadronic 
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decays and the latter two 6-tags within a fat jet or two separate slim ones. Some of these 
channels were also studied in the comprehensive review of exotic Higgs decays contained 
in |58]. 

In the present paper, working under the assumption that a light Ai state has been found 
through one or more of the decay channels analysed in |49) . we assess the scope of the LHC 
Run H in profiling its nature by resorting instead to the VBF and HS Higgs production 
modes. This would then enable access to the heavy Higgs couplings to both charged and 
neutral gauge bosons, thus complementing the GF channel which only allows one to measure 
their fermion couplings. It is worth emphasising that for the non-SM-like scalars, this might 
be the only chance to access these couplings as the decay to pairs of pseudoscalars can be 
completely dominating. Furthermore, although the VBF and HS channels have significantly 
smaller cross sections than GF, the improved possibilities for tagging through the additional 
(forward/backward) jets for VBF and the additional leptons from vector boson decays in 
HS, may still render them competitive against GF, for which many of the triggers needed 
for fully hadronic final states have been scarcely tested in the experimental environment. In 
view of this, in this paper we perform a preliminary study of light His in the VBF and HS 
Higgs production modes by making some simplifying assumptions, in particular by ignoring 
detector effects and trigger thresholds. 

Before plunging into the details of this new analysis, we should also point out here that 
in the NMSSM both Hi and H 2 , the lightest and next-to-lightest GP-even Higgs bosons, 
respectively, can alternatively play the role of the SM-like Higgs boson PfsM) as emphasised 
already in |59l [ 6 OI El EJ and confirmed in jlHl |56l • 

The article is organised as follows. In section we define the parameter space of the 
NMSSM under consideration. In section we explain our approach to scan the NMSSM 
parameter space while in section]^ we describe our signal-to-background analyses. Then in 
sections and we discuss in detail our results for VBF and HS (the latter separately 
for the neutral {ZH) and charged (WH) channels) at the LHG over the entire NMSSM 
parameter space. Then, after presenting some benchmark points available for experimental 
investigation in section and testing them at the GERN pp machine, we afford a brief 
digression on the physics of a light Ai state at a future e'’'e“ collider (in section]^. We 
summarise and conclude in section [101 

2 The NMSSM parameter space 

The idea behind the NMSSM is to explain the peculiar feature of the MSSM that the 
supersymmetry preserving p term is phenomenologically required to be of the same scale as 
the soft supersymmetry breaking mass parameters, while in principle it would be expected 
to be of a completely different origin. 

In the NMSSM this so-called p problem is solved by introducing an extra gauge-singlet 
chiral superfield S whose scalar component receives a vacuum expectation value (VEV) 
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due to its soft supersymmetry breaking terms. All that is needed to generate an effective 
/i term of the same order as the soft supersymmetry-breaking parameters, is then to have 
a term XSHuHd in the superpotential and the (effective) /i term = As will be given by 
the VEV of S times the coupling constant A. We also need to add a cubic term in S to the 
superpotential, so that the terms involving S read 

XSHuHd+'^S^ ( 1 ) 

where A and k are dimensionless coupling constants. Further, one needs to add the corre¬ 
sponding soft supersymmetry breaking terms in the scalar potential. The soft supersym¬ 
metry breaking terms relevant for the Higgs sector are: 

+ i^^xSHuHd + —AkS^ -\- h.c.), ( 2 ) 

where and are dimensionful mass and trilinear parameters. By 

minimising the scalar potential we can trade the three scalar mass parameters for |Ueff and 
tan/3 (i.e., the ratio between the VEVs in the up type and down type Higgs doublets, 
Vu/vd)- We also replace with Mp, the singlet pseudoscalar mass entry; if input is given 
at the Electroweak scale, this is close to the lightest pseudoscalar mass which improves the 
efficiency of the scan. For this reason we use parameters dehned at the EW scale throughout 
the paper. 

For the sfermion masses we use a common mass parameter mo and for all gaugino 
masses we use a common parameter m ^/2 but to take into account the effects of running 
from a high scale unihcation we use M 2 = rnij 2 , Mi = mi/ 2 / 2 , M 3 = 3 mi/ 2 - Similarly, we 
use a common parameter Aq for all trilinear parameters except A;^ and A^. Since the most 
signihcant effect these parameters have on the Higgs sector is radiative corrections to the 
Higgs mass, unifying them in the above manner should not impact much on the observables 
of interest, thus giving maximum freedom in the Higgs sector while keeping the number of 
free parameters at a manageable level. This leaves us with nine parameters: 

mo, mi/ 2 , Ao, tan/3, /ieff, A, k. Ax, Mp. (3) 

As can be seen in the tree-level mass for the lightest doublet scalar, i.e. essentially (up 
to mixing with the singlet scalar) UgM in the NMSSM |T], 

~ m| cos^ 2/3 -|- A^u^ sin^ 2/3 , (4) 

there is an additional contribution to the Higgs mass coming from the A term, not present 
in the MSSM. This means that, when this term is sizeable, i.e., when A is large and tan/3 
is small, one can obtain the measured 125 GeV mass without resorting to large radiative 
corrections from the stop sector, which is necessary in the MSSM. In the forthcoming 
sections we will refer to this part of parameter space as the ‘naturalness limit’. Also notice 
that, when H 2 = ^^SM, it is possible that some mixing between Hi and H 2 increases 
mH 2 further, hence making it even easier to achieve (and indeed exceed) the required 125 

GeV [Sni EH E3 [Ml EH EH]- 
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3 The scans 


In order to investigate the prospects of the discnssed channels in the NMSSM we per¬ 
formed a nnmber of Bayesian scans of the above mentioned 9-dimensional parameter space. 
These scans nse MnltiNest-v2.18 [69] for nested sampling of the parameter space and 
NMSSMTools-v4.6.0 |70], inclnding the bng-fix of version 4.7.1, for the calculation of the 
Higgs mass spectrum, couplings. Branching Ratios (BRs) and constraints on the parame¬ 
ter points (including LEP-II bounds on light scalars and pseudoscalars). The output from 
NMSSMTools is further processed by HiggsBounds-v4.1.3 imiTIllTallTl] to ensure against 
exclusion from searches for other Higgs bosons. Also SuperIso-v3.3 m is used to calculate 
6 -physics variables. These are then required to comply with the constraints from m (in 
all cases, the last number refers to a theoretical uncertainty in the numerical evaluation): 

• BR (Bs ^ n+fj,-) = (3.2 ± 1.35 ± 0.32) x 10“^ 

• BR (Bu ^ Tiy) = (1.66 ± 0.66 ± 0.38) x 10"^, 

• BR (B ^ Xs 7 ) = (3.43 ± 0.22 ± 0.21) x lO”'^. 

To guard against over-closure of the Universe, an upper bound of < 0.131 on the 
dark matter relic density was also applied. This was set assuming a 10% theoretical error 
on the central value of 0.119 from PLANCK m and the relic density was calculated with 
the help of MicrOMEGAs-v2.4.5 |78j . We do not include any constraint on {g — 2)^; the 
contributions here are always too small to cause problems, they are actually too small to 
reconcile the experimental value with the Standard Model calculations. 

A total of four scans were run; two where either one of Hi and H 2 was allowed to be 
the SM-like one, and since this mostly gave points with Hi = PfsMj two other scans were 
run including a bias towards H 2 = H^m- In both cases did we use one scan focusing on the 
naturalness limit and one wider scan. The parameter ranges for all scans are given in table[^ 
where the reduced range focuses on the naturalness limit and, since the couplings relevant 
to e.g.. Hi —7- AiAi decays depend on A this is the most promising part of parameter space 
to look for these channels. All the scan results are then combined into two samples, one 
with Hi = Hsm and one with H 2 = PfsM- 

For all samples we require to lie between 122 and 129 GeV. The best experimental 

values of the Higgs mass are 125.03 GeV from GMS [79] and 125.36 GeV from ATLAS [80] 
with uncertainties of the scale of a fraction of a GeV, however, to allow for potentially 
large theoretical uncertainties, we allow a significantly larger range. Note, though, that the 
benchmark points of section]^ are within the experimental limits. 

The scans also contain a bias towards low pseudoscalar masses (favouring ^ 65 
GeV, but allowing up to 140 GeV) and SM-like signals rates for PfsM- These con¬ 
straints are implemented to ensure that the scans do not waste too much time exploring 
uninteresting parts of parameter space. Gonversely, they should not exclude any points 
that might be of interest for further investigation. 


6 


Parameter 

Extended range 

Reduced range 

mo (GeV) 

200 - 2000 

200 - 2000 

mi /2 (GeV) 

100 - 2000 

100 - 1000 

Alo (GeV) 

-5000 - 5000 

-3000 - 3000 

/ieff (GeV) 

50 - 1000 

100 - 200 

tan (3 

1 - 30 

1 - 6 

A 

0.01 - 0.7 

0.4 - 0.7 

K 

0.01 - 0.7 

0.01 - 0.7 

AIa (GeV) 

200 - 2000 

200 - 1000 

Mp (GeV) 

3 - 140 

3 - 140 


Table 1: Parameter ranges used in the scans. The reduced range focuses on the naturalness 
limit. 


Though similar scans were also performed in our previous paper [19], there are a num¬ 
ber of important updates. Apart from updated versions of some of the software (most 
importantly NMSSMTools), the most notable difference is the treatment of the constraints 
on signal rates for Hsm] this is now done using the built in constraints of NMSSMTools, 
which uses constraints obtained from Lilith 1.0 |81) . Especially important is the inclusion 
of signal constraints on H^m bb] it turns out that virtually all points in the naturalness 
limit with niAi < excluded by this constraint. Given the poor measurement of 

this channel as compared to e.g. iLsM 77 and PfsM ZZ, this might sound surprising, 
but it turns out that the latter channels can afford a relatively large BR(LfgM —?■ AiAi) by 
compensating with larger branching ratios for the measured channels, this however is not 
possible for /Lsm bb since this branching ratio cannot be further increased and points 
with increased vector boson channels tend to have some singlet component in the /LgM 
which further reduces the production cross-section and therefore conflicts with H^m bb. 
In the case of Hi = LfgM this only conhrms the picture from where only a wider scan 
yields points below H^yi/2. For H 2 = on the other hand, this makes a big difference; 
while our previous studies deemed Hi —)• AiAi with H 2 = PfsM a very promising channel, 
this has now changed. 

To illustrate the difference, we plot in hgure the same sensitivity curves as in hg- 
ure 12(a) of |49| but with our new scans instead. It is clear that, while the earlier scan 
found good chances of discovery already at 30/fb the prospects are now looking far dim¬ 
mer. If we look at A as a function of niAi as is done in hgure we see clearly that it is 
the high A region that is unreachable for mA^ < in this region the precence of 

the i^sM channel suppresses the signal rates of H^yi, especially for large A (the 

relevant H^mAiAi coupling contains a factor A^); for the vector boson channels this can 
be somewhat compensated by giving a somewhat larger singlet component to the i^sMi 
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Figure 1: The sensitivity in the previously most promising channel with gg —)• Hi —)• 
with H 2 = F^sM) sensitivity curves are identical to the ones in figure 12(a) of |49| . The 
sensitivity curves correspond to the one fat jet plus two r-jets (the curves at low ruAi ) and 
the two single b-jets plus two r-jets (the curves at higher niAi) analyses. For the fat jet 
plus two r-jets rriHi = 100 GeV, while for the two b-jets plus two r-jets curve rriHi = 125 
GeV, which allows the coverage of points with large • The green points comply with all 
constraints of our scans. The red triangles mark the benchmark points defined in section Isl 


this reduces the production cross section but also the width of competing decay channels 
hence increasing the signal rate for the vector boson channels; however, this cannot be done 
simultaneously for FFsM bb. With these updated constraints BR(FFsm ^ 1 ^ 1 ) never 
reaches more than just above 20 %. 

It is also worth mentioning that compared to the scans of |49) . we here use somewhat 
different parameters, especially, we use the singlet pseudoscalar mass term, Mp, as input 
(to improve efficiency) and all input is at the EW scale (a necessity if Mp is to actually 
correspond to mA^)- We also use the full NNLO contributions to the scalar masses in 
NMSSMTools, in contrast to what was done in |49| . This raises the especially for 

large stop mixing and mass, and might therefore further favour points with large mnsM 
contributions from the stop sector rather than the NMSSM specific contributions of the 
naturalness limit. 

Finally, in both figure |^a) and (b), there is a depletion of points above 70-80 GeV: this 
is due to the bias towards low niAi included in the scans. 

Given the requirement of SM-like signal rates for Hsm, it is natural that its reduced 
coupling to vector bosons will always be large (close to 1). The reduced coupling of Hs will 
be smaller but can, due to mixing, still be significant, however the reduced coupling of FF 3 
always turns out very small. 

To understand this let us look at the couplings of interest: they both have a factor [1] 

{VdS31+VuS32), (5) 
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Figure 2: The parameters A versus ruAi- The colour code for the points is the same as in 
figure 


where S 31 and S 32 are elements of the neutral scalar mixing matrix, defined as Hi = 
with We know that Vu/vd = tan/3 by definition 

and, due to the structure of the mixing matrices, it turns out that S' 3 i/S '32 ~ — tan/3. 
Hence the factor (|^ becomes small; typically, the reduced coupling tends to be ~ 0.01 
or smaller and, since it comes squared in the production cross sections for both VBF and 
HS, we get a suppression factor of at least 10“^, rendering these channels useless for H 3 . 
Note that the cancellation mentioned above gives a much stronger suppression than naively 
expected from the fact that H^m has to have SM-like couplings combined with the sum-rule 
requiring the sum of the squares of the vector boson couplings for all three scalars to be 
unity. 

4 Event analyses 

Since the vector boson couplings to H 3 are always very small we do not consider channels 
including H 3 production here, they do not have big enough cross sections to be of any 
interest. Furthermore, since the channels Hi ^2 were shown in |49| to be very 

difficult, we here focus only on the Hi ^2 channels. Given the lower cross sections 

for VBF and HS production as compared to GF, none of the other channels carry promise 
for detection. 

To estimate the sensitivity in the channels of interest we use MadGraph5_aMG@NLO |82) . 
employing default parton distribution functions and factorisation and renormalisation scales, 
to generate the relevant backgrounds. Hadronisation and signal generation is then done 
using Pythia 8.180 [83], while jet clustering and jet substructure studies are done using 
FastJet-v3.0.6 |84j . 
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The production cross sections for the signals are calculated using tabulated cross sections 
for the SM Higgs, taken from |85| . together with reduced couplings for the relevant scalar 
from NMSSMTools. 

For all final state objects we use the following acceptance cuts: 

• \r]\ < 2.5 for all final state objects, 

• px > 15GeV for all final state jets (r-jets, 6-jets and light-quark jets), 

• Pt > lOGeV for all final state leptons (e^, 

• Ai? = yj (A? 7)2 -|- (A(/))2 > 0.2 for all 6-quark pairs, 

• Ai? > 0.4 for all other pairs of final state objects, 

where pT, rj, (j) are the transverse momentum, pseudorapidity and azimuthal angle, respec¬ 
tively. Both H-mesons and rs are decayed by Pythia and the decay products included in 
the jet clustering in Fast Jet. The identification of 6- and r-jets is then done by tracing the 
jet constituents back in the event record, thereby identifying 6- and r-jets as jets where 
most of the constituents are decay products from H-mesons and rs respectively. On top of 
this the event sample is rescaled by an assumed tagging efficiency of 50% for each 6- and 
r-jet. This efficiency is somewhat conservative (c.f. |86[ I87j l but this is motivated by our 
soft signals. No knowledge about the charge of the jets is assumed. For the ti background 
for the WH channel, we assume a 1% mistagging probability for light jets. 

For the VBF backgrounds the two additional light jets are required to satisfy the fol¬ 
lowing criteria: 

• H < 5 , 

• Pt > ^0 GeV, 

• AR>3, 

• Mjj > 300 GeV, 

where Mjj is the invariant mass of the two forward/backward jets. (Due to numerical 
difficulties in producing enough statistics, somewhat harder cuts, Mjj > 500 GeV and 
Pt > 40 GeV, are imposed on the parton level production in the 46 final state.) In addition, 
the jets are required to be located in the opposite hemispheres (i.e. one positive and one 
negative r]) and there should be no other jet (apart from the signal objects) with pT > 30 
GeV between (in terms of r]) the two jets. 

The cross sections for the dominant backgrounds after the above acceptance cuts are 
applied, are given in table Note that for the Oj and Ij backgrounds the above VBF 
specific cuts are not applied, hence the much larger cross-sections, this will be compensated 
by smaller acceptance after hadronisation. After hadronisation and jet clustering the final 
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Ghannel 

Parton level cross section 

VBF 

2j + 46 

72 pb 

Oj -|- 262r 

3.1 pb 

lj + 2b2T 

3.2 pb 

2j + 262r 

0.19 pb 

Oj -|- tt 

597 pb 

Ij + ti 

845 pb 

2j + tt 

80 pb 

ZH 

Z + 46 

0.31 pb 

Z -F 262t 

2.7 fb 

WH 

IT+ 46 

36 fb 

ttbb 

4.0 pb 

tt 

597 pb 


Table 2: Background cross sections for the dominant backgrounds at parton level as calcu¬ 
lated by MadGraph. 


VBF cuts are applied to the real jets (as opposed to the parton as done above) and then 
Mjj > 500 GeV and px > 40 GeV are required for all signal and background in the VBF 
channels. For the rest of the paper, in all discussions of VBF production tt refers to the 
sum of Oj -|- tt, Ij + tt and 2j + tt and 262r refers to the sum of Oj -|- 262r, Ij + 262t and 
2j + 262t. The total background for the 262r channel is the sum of tt and 262t. We do 
not include further backgrounds in the 46 channel since this channel is in any case inferior 
to 262r and therefore of little interest. 

In order to optimise the sensitivity to boosted Ais we employ the jet substructure 
method of [88] (see also m for further details). This gives us fat jets that we assume to 
originate from an Ai decaying into a 66 pair. To avoid contamination from single 6-jets, we 
require all fat jets to have pT > 30 GeV and invariant mass > 12 GeV. 

For each event with the proper hnal state we look for two Ai candidates (either one fat 
jet, two 6-jets or two r-jets depending on which channel we are looking at) and compare their 
respective invariant mass: if they are within 15 GeV of each other we accept that event and 
if, in addition, the combined invariant mass of the two candidates is within 125 ± 30 GeV, 
it is accepted as an event where an H^m was produced and decayed into two ^is. In all the 
channels we perform two analyses in parallel: one with the jet substructure method, where 
only fat jets and no single 6-jets are used, and one where no jet substructure is exploited 
but only single 6-jets are used. For the four 6-jet hnal state we check all combinations of 
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6-jet pairs and accept the first one with both invariant masses within 15 GeV of each other. 

To obtain the sensitivity for a given rriAi we then count the events where the Ai can¬ 
didates masses are within 15 GeV of the mass mA^ and can then calculate — for a given 
luminosity, £ — how large a signal cross section is needed to obtain a significance S/y/B > 5 
where S is the number of signal events and B is the background. In all channels we require 
at least 10 events in order to claim discovery, so if S/y/B > 5 is fulfilled for S < 10, we 
instead use 5* > 10 as the limiting sensitivity. This means that, in channels with very low 
background, the sensitivity is oc £ rather than oc y/C as is the case for S/y/B. 

We will be studying the VBF and HS channels in turn, splitting the latter into the ZH 
and WH modes. For the sensitivity curves in the forthcoming sections, we limit ourselves 
to scalar masses up to 175 GeV, there is no problem in principle to employ these analyses 
at higher masses, they would actually gain some sensitivity with harder objects to study, 
but the production cross-section for the initial scalar also drops fast with increasing mass, 
hence prohibiting any discovery prospects. 

5 VBF 

The backgrounds used in the analysis of the VBF mode are irreducible, i.e., 46-1-2 jets and 
262r-|-2 jets. In addition, we include the tt+2 jets — with both VFs from the top quarks 
decaying to rs — in the 262r channel. Since the latter turns out to be the dominant 
background in this channel, we also invoke a cut on Missing Transverse Momentum (MET) 
to reduce it. 

This cut requires that the px of the two r-jets combined should be larger than the total 
MET of the event. This reduces the tt+2 jets background by a factor 2-3 while leaving 
both signal and irreducible backgrounds virtually intact. Note, however, that the MET 
here is simply the sum of the momentum of all invisible particles (i.e., neutrinos), a full 
detector simulation with mis-measured/missing jets, pile-up etc. would be necessary to fully 
determine the true effectiveness of this cut. This latter cut, though effective, is, however, 
not crucial for the usefulness of this analysis. 

Figure]^ shows the discovery reach in the interesting channels for 3000/fb of integrated 
luminosity and using the overall constraint that the total four-body invariant mass should 
be 125±30 GeV. It is clear that the 262r channels are the most promising ones and hence we 
will not consider the 46 channels in the following. All sensitivity curves have been rescaled by 
a factor 1/0.9 for each bb pair in the final state and a factor 1/0.1 for each rr pair in the final 
state, to allow for a direct comparison of the sensitivity to cr{qq —>■ qqHi) x BR(i7j —)• AiAi). 
As expected, we see in figure that the jet substructure analysis only works well for rather 
low ruAi ■ 
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Figure 3: Sensitivity for the 46 and 262r channels in VBF production. All curves assume 
3000/fb of integrated luminosity and use m//. = 125 GeV and the corresponding cut on the 
total final state invariant mass. 

5.1 Hi SM-like 

To illustrate the reach of the LHC in these channels, in figure]^ we compare the sensitivity 
with the points from our scans when Hi = Hqm- In figure |^a), we show sensitivity curves 
using a 125 GeV Higgs (and employing the corresponding constraint), always using the best 
out of the 2 6-jets and fat jet analyses. As is clearly seen, the prospects for Hi —)• AiAi 
will always be rather limited due to the small number of points in the interesting region of 
parameter space. It is clear that the points never reach above 1 pb, this is a direct effect 
of the constraints on the signal rates of H^m, higher BR(L6gM — AiAi) would necessarily 
suppress other channels unacceptably much. This means that we need 3000/fb to discover 
something and our best chance would be the low mass region where the jet substructure 
methods can improve sensitivity. 

As regards the H 2 —)• AiAi channel, it is somewhat easier to find acceptable points just 
above the ttia^ < threshold, however, they do not reach much higher in terms of 

rates. This can be seen in figure |^b), where the prospects for H 2 —)• AiAi are illustrated. 
Figure |^b) uses sensitivity curves for = 175 GeV in order to cover the whole inter¬ 
esting parameter space. Also, here do we always employ the analysis (with or without jet 
substructure) with the best sensitivity. The smallness of the rates also above the kinematic 
threshold is a consequence of the requirement that TLsM has to have very SM-like couplings 
to meet the signal rate constraints, and hence H 2 has to be very singlet-like and hence have 
very small production rates even if the Hi —)• AiAi channel is kinematically closed. In the 
end this means that our best chance also in this channel seems to be the light pseudoscalars 
where jet substucture is useful. 
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Figure 4: LHC reach in Hi —>■ AiAi (left) and H 2 —>■ (right) for Hi = H^yi in 

the VBF channel. In panel (a) the sensitivity curves use =125 GeV, while (b) uses 

mH 2 = 175 GeV and both uses the best of the 2 6-jets and fat jet analyses. The colour code 
for the points is the same as in hgurej^ 


5.2 H 2 SM-like 

As mentioned before, after the inclusion of especially the signal rate constraint on TLsM 
66, the case H 2 = Hsm does not differ so much from Hi = Hsm- In hgurej^a) we show the 
sensitivity for Hi —?• AiAi, both with jet substructure using vriHx = 100 GeV and without 
jet substructure using tuhi — 125 GeV, but without constraining the four-body invariant 
mass. We use two separate curves since they use different The use of = 100 

GeV is motivated by the typical values of while 125 GeV is needed to cover the whole 
range of niAx , the = 100 GeV curve would be kinematically cut off at = 50 GeV. 
Also in this channel we will need 3000/fb to see something. 

From hgure [^a) it even looks like the jet substructure does worse than the simpler 2 
6-jet analysis, though one should remember that these curves uses different the lower 

ruH^ of the fat jet analysis means decreased sensitivity. For the 46 hnal state this can be 
compensated by an improved sensitivity due to the fat jet analysis, however, as can be seen 
in hgure the low mass improvement with the fat jet analysis is not as great for 262 t as 
it is for 46, leading to comparatively poor sensitivity to ~ 20 GeV in hgure j^a) . 

In hgure j^b) we show the sensitivity in the H 2 —)• AiAi channel using mH 2 = 125 GeV 
and requiring the four-body invariant mass to be 125 ± 30 GeV. For all points along the 
sensitivity curves we use the analysis that gives the best sensitivity (with or without jet 
substructure). As expected, this is very similar to hgure |^a) with a clear upper limit for 
the rate at 1 pb stemming from the signal rate requirements on iLsM- 
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Figure 5: LHC reach in Hi —)• AiAi (left) and H 2 —)• (right) for H 2 = Hsm hr the 

VBF channel. The sensitivity curves in panel (a) correspond to the one fat jet (the curves 
at low mAi using rriHi ~ GeV) and the two single b-jets (the curves at higher niAi 
using rriH^ = 125 GeV) analyses. In panel (b) mH 2 — 125 GeV is used along with the best 
of the 2 6-jets and fat jet analyses. The colour code for the points is the same as in hgurej^ 

6 Higgs-strahlung via ZH 

With an additional Z boson, triggering and background suppression is much improved. To 
extract the signal we only use leptonically decaying Z bosons. This means, that in addition 
to acceptances, we require one di-lepton pair with invariant mass 90 ± 10 GeV. This is very 
powerful in suppressing the backgrounds, however, the small leptonic BR of the Z together 
with the small production cross sections in the ZH channel, means that one will struggle 
to get a large enough signal. As one would then expect, the best hnal state to look for 
is not 262r as was used before, but rather 46 which gives the highest signal rate. This is 
clearly seen in hgurej^ where the sensitivity in the various channels are shown for 3000/fb 
of integrated luminosity. All channels use m//- = 125 GeV and the corresponding cut on 
the total hnal state invariant mass. 

6.1 Hi SM-like 

Similarly to the VBF case, the SM-like Hi scenario is difficult with respect to detection. 
In hgure [^a) we show the sensitivity in the Hi —)■ AiAi channel, given ruHi = 125 
GeV (with a corresponding cut on the four-body invariant mass), as compared to the 
acceptable parameter points. We notice that the prospects are even dimmer than in the 
VBF case, basically our only hope is the jet substructure methods that gives a rather 
large improvement for low for the 46 hnal state used here this improvement can be 

signihcantly larger than for the 262r used in the VBF channels. As usual only the analysis 
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Figure 6: Sensitivity for HS off a Z boson. All curves assume 3000/fb of integrated lumi¬ 
nosity and use rriff. = 125 GeV and the corresponding cut on the total final state invariant 
mass. 



Figure 7: LHC reach in Hi —)• AiAi (left) and H 2 —>• AiAi (right) for Hi = TLsm in the ZH 
channel. In panel (a) the sensitivity curves use rriHi =125 GeV, while (b) uses = 175 
GeV and both uses the best of the 4 6-jets and fat jets analyses. The colour code for the 
points is the same as in figure 


with the best sensitivity is used in each point of the curves. 

Moving on to H 2 —)• AiAi, we see in figure j^b) that there is only marginal hope for 
detection even at 3000/fb. The curves here correspond to mH 2 = 175 GeV, with no cut on 
the total final state invariant mass and always using the most efficient of the single 6-jets/fat 
jets analyses. The main reason for the difficulties in this channel is clearly the small signal 
rates, but also the absence of a cut on the overall invariant mass impacts negatively on the 
sensitivity. 
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Figure 8: LHC reach in Hi —)• AiAi (left) and H 2 —)• (right) for H 2 = Hsm hr the 

ZH channel. The sensitivity curves in panel (a) correspond to the two fat jets (the curves 
at low using niHi ~ GeV) and the four single b-jets (the curves at higher niAi 
using ruHi — 125 GeV) analyses. In panel (b) = 125 GeV is used along with the best of 

the 4 5-jets and fat jets analyses. The colour code for the points is the same as in figure 

6.2 H 2 SM-like 

With the lower sensitivity as compared to VBF production, the H 2 = L^SM is also rather 
difficult. As can be seen for Hi —?■ AiAi in figure |^a) — where the sensitivity curves use 
ruHi — 100 GeV with jet substructure, and ttihi — 1^5 GeV with single 5-jets (no constraint 
on overall invariant mass), respectively — even 3000/fb is probably not sufficient for any 
discovery. 

In figure j^b) we show the reach for H 2 —t- AiAi, the sensitivity curves correspond to 
mH 2 = 125 GeV and the corresponding constrain on the four-body invariant mass, always 
showing the most sensitive analysis. As before, this is basically identical to figure j^a), with 
only marginal discoverability for low masses where jet substructure may be useful. 

7 Higgs-strahlung via WH 

In general WH production will have higher cross sections than ZH production and, since 
the W also has roughly three times higher leptonic BR as compared to the Z, this channel 
will exhibit much higher rates. In order to tag the W we require exactly one isolated lepton 
in the event, in addition to the signal objects. Similarly to ZH production, we only look 
at the channel with the highest rate, i.e., W + 45. 

There has been a number of earlier studies of the VF -|- 45 channel, including parton 
level studies |89| and |90| and a full detector study in |91) . While the parton level analysis 
of [90] arrived at significantly higher sensitivity than we did, our results are in reasonable 
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agreement with |91) . (Note however that, although both of these studies use four 6-tags as 
well as a cut on the four-body invariant mass, neither of them uses the requirement that 
both Ai candidates should have similar mass.) 

In addition to the irreducible -|- 46 background, there are signihcant backgrounds 
from tt — with two light jets from one of the resulting Ws being mistagged as 6-jets — 
as well as from ttbb events. The latter is in our studies the most signihcant one, often at 
least one order of magnitude larger than the irreducible background. This conclusion is in 
agreement with |90) . while |91) hnds that detector smearing pushes the tt background to 
lower invariant masses and hence becomes a signihcant background also at 125 GeV. To 
suppress the ttbb background we employ a veto against hadronically decaying Ws. This 
means an event with two light jets with pT > 15 GeV and combined invariant mass = 80±15 
GeV is rejected as it is likely to come from a ttbb event with one W decaying leptonically 
and the other one hadronically. The events where both Ws decay leptonically should be 
suppressed by the fact that we ask for exactly one lepton and hence reject events with two 
isolated leptons (this also suppresses any Z + 4:b backgrounds). 

Since the smallness of the signal rates in both WH and (even more) in ZH production 
is a bigger problem than background suppression, one could consider requiring only three 
6-tags. This was the approach of j89j and should yield signihcantly higher signal rates, 
especially as one 6-jet is often missed due to pT cuts, etc. However, this requires a much 
more detailed study as there are many more contributing sources of background. It is 
also not clear how to implement the invariant mass constraints as one have to assume 
that sometimes the fourth 6-jet is not selected even as a light jet. Such considerations are 
therefore beyond the scope of this paper. 

7.1 Hi SM-like 

In this case higher rates as compared to ZH production does mean better discovery 
prospects even though the background is also larger due to ttbb. 

In hgure a) we show the discovery reach in the Hi —)• AiAi channel in the Hi = UsM 
case. The sensitivity curves are set for tuhi = 125 GeV and the corresponding cut on the 
overall invariant mass. If we compare with hgure [^a), the reach here is much greater, even 
300/fb might be enough for detecting some His around 20 GeV. Such a discovery is not 
possible even for VBF production, as can be seen in hgure j^a), where 300/fb does not reach 
the upper 1 pb limit for the rate. The reason for the relative success oiWH production for 
these low masses is that in the 46 hnal state the use of fat jets leads to a more signihcant 
improvement over a single 6-jet analysis than is the case for the 262r channel used for VBF 
production. Also, the backgrounds in the WH channels are more severe at higher masses, 
rendering the low mass region more interesting. 

The H 2 —7- AiAi channel is much less optimistic, as can be seen in hgure j^b), where 
we set mH 2 = 175 GeV with no cut on overall invariant mass for the sensitivity curves. It is 
clear the growth of the ttbb and ti backgrounds with increasing invariant masses suppresses 
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Figure 9: LHC reach in Hi —>■ AjAi (left) and H 2 —>■ (right) for Hi = Hqm hr 

the WH channel. In panel (a) the sensitivity curves use mn^ =125 GeV, while (b) uses 
mH 2 = 175 GeV and both uses the best of the 4 6-jets and fat jets analyses. The colour 
code for the points is the same as in hgure 


the sensitivity at higher masses. One could try some invariant mass cuts to suppress these 
backgrounds and possibly reach the points just above the mA^ = threshold, but we 

leave such considerations to section IH 

7.2 H 2 SM-like 

The case H 2 = is rather similar to the Hi = L7sm case. As can be seen in hgure [To]) a) 
— displaying sensitivity curves with = 100 GeV (using fat jets) and 125 GeV (using 
6-jets) but with no overall cut on invariant mass — there is just a marginal hope for 
detection in Hi —)• AiAi at 3000/fb. We do see improved sensitivity in the low mass 
region as compared to VBF production, but not enough to cover any signihcant part of the 
parameter space. This is again a consequence of the strength of the fat jet analysis for 46 
as well as the main backgrounds tibb and ti being smaller for lower invariant masses. 

For H 2 —)• AiAi (hgure [T^b)) we again see essentially the same as hgure [^a). Here 
we use mH 2 = 125 GeV and constrain the four-body invariant mass of the hnal state to be 
125 ± 30 GeV. The sensitivity curves use a combination of the fat jet and the single 6-jet 
analyses, always showing the more sensitive one. 
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Figure 10 : LHC reach in Hi —)• AiAi (left) and H 2 —)• (right) for H 2 = H^m hr 

the WH channel. The sensitivity curves in panel (a) correspond to the two fat jets (the 
curves at low rriAi using rriH^ = 100 GeV) and the four single b-jets (the curves at higher 
rriA^ using = 125 GeV) analyses. In panel (b) = 125 GeV is used along with the 

best of the 4 6-jets and fat jets analyses. The colour code for the points is the same as in 
hgure 


8 Benchmark points 

As mentioned before, the detection sensitivities in all the channels studied in this paper 
are signihcantly worse than the corresponding results from GF productionj^ Furthermore, 
this is generally true for all parameter space regions that we have been able to access in 
the analyses carried out here. Therefore it is fair to assume that, by the time the channels 
under investigation become interesting for experimental study, the lightest scalar {Hi and 
H 2 ) as well as pseudoscalar (Ai) Higgs states will already have been discovered via GF and 
the goal for both the VBF and HS channels will become to enable one further study their 
properties, as explained in the introduction. 

For this purpose we define three benchmark points and study them under the assumption 
that we know the masses involved and can hence use this information to further constrain 
the kinematics in the attempt to increase the sensitivity]^ The details of the three points 
are given in table The points are chosen to cover as much as possible of the parameter 
space within reach. Point 1 has a rather light Ai of 21 GeV and hence lies in the region 
where jet substructure methods are of importance, while point 2 is closer to the threshold 
for Hi ^2 AiAi with = 55 GeV. The last point is the only one with Hi = /Lsm 

®This is true unless the triggering in the GF channel turns out to be more challenging than what is 
presently hoped for. 

'*Even without prior knowledge of the masses, these analyses can be employed by scanning over the 
masses involved. This, however, would render large look-elsewhere effects. 
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Gase 


= H2 

-f^SM = Hi 


Point 1 

Point 2 

Point 3 

Input parameters 


niQ (GeV) 

1516.0 

1470.5 

1326.7 

mi /2 (GeV) 

491.58 

947.19 

369.23 

^0 (GeV) 

-2966.1 

2523.0 

2839.5 

f^es (GeV) 

113.54 

122.98 

189.28 

tan /3 

13.051 

12.46 

2.882 

A 

0.08899 

0.1641 

0.5408 

K 

0.04042 

0.07769 

0.1818 

AIa (GeV) 

781.94 

1714.9 

479.29 

Mp,(GeV) 

21.77 

59.38 

68.52 

Observables 


ruAi (GeV) 

21.12 

54.87 

63.26 

(GeV) 

100.54 

110.77 

125.45 

rriH^ (GeV) 

125.22 

125.31 

138.90 


Table 3: Some specifics of the three benchmark points. 


and with = 63 GeV, it is designed for H 2 —)• AiAi studies. Below the —)• 

threshold, the phenomenology of Hi = H^yi and H 2 = H^m is very similar, obviating any 
need for additional benchmark points in that region. 

in the study of these benchmark points we therefore use a somewhat modified kinematic 
analysis. Specifically, all Ai candidates (i.e., 6-jet pairs, fat jets or r-jet pairs) are required to 
be within 15 GeV off the (assumed known) Ai mass. For each point we run two simulations, 
one for Hi —?■ and one for H 2 —?■ AiAi and in each case the combined invariant mass 

of the two Ai candidates is required to be within 30 GeV of the (assumed known) Hi or 
H 2 mass. 

The result of these studies are displayed in table where the cross sections after all 
cuts are presented for the signal as well as the backgrounds. We also show the integrated 
luminosity needed to obtain S/y/~B > 5 with at least 10 events. The result of the jet 
substructure (fat jet) analysis is only shown for point 1 as this is the only scenario with 
an Ai light enough for such studies to be useful, though, in that case, this is clearly the 
most effective approach. Also, we do not include Hi —?• AiAi for point 3 as this channel is 
kinematically closed. 

As stated before, the signals are in general larger than the backgrounds and, in many 
cases, the main constraint on the required luminosity is the requirement of at least 10 
events. We also see in table that WH \s usually the most promising channel as it has 
the highest rates. It also has the highest backgrounds though, due to tibb and ti and, since 
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Point 1 

Point 2 

Point 3 


H = Hi 

II 

II 

tq“ 

H = H2 

II 

tq“ 


6-jets 

fat jet 

6-jets 

fat jet 

6-jets 

6-jets 

6-jets 

qqHi —!■ qq2b2T 


Signal [pb] 

2j + 262r 

2j + tt 
£ [fb-i] 

1.7 X 10-*’ 
2.4 X 10-'^ 

6.8 X IQ-'^ 

8400 

2.7 X 10-5 
8.5 X 10-'^ 
2.2 X 10-5 
10000 

4.1 X 10-5 
2.4 X 10-'^ 

6.8 X 10-'^ 

2400 

8.8 X 10-5 

8.5 X 10-'^ 
2.2 X 10-'^ 
1100 

3.9 X 10-5 
2.6 X 10-5 
1.1 X 10-5 
5800 

7.0 X 10-5 
2.6 X 10-5 
1.1 X 10-5 
1900 

1.0 X 10-5 

3.4 X 10-5 

1.5 X 10-5 
1200 

ZH, 2^46 


Signal [pb] 

Z + ib [pb] 

£ ]fb-i] 

1.4 X 10“® 
9.1 X 10“’^ 
11000 

4.1 X 10-5 

1.1 X 10-5 
2400 

2.9 X 10-5 
8.7 X 10-'^ 
3500 

7.8 X 10-5 
1.1 X 10-5 
1300 

2.2 X 10-5 
2.7 X 10-5 
13000 

2.9 X 10-5 

3.9 X 10-5 
11000 

4.0 X 10-5 
4.0 X 10-5 
6300 

WHi £46 


Signal [pb] 

IT -f 46 ]pb] 
ttbb [pb] 
ti [pb] 

£ ]fb-i] 

1.1 X 10-5 

1.1 X 10-5 

8.2 X 10-5 

9.3 X 10-'^ 
2300 

3.9 X 10-5 
1.7 X 10-5 
1.4 X 10-5 

4.9 X 10-'^ 
266 

2.1 X 10-5 

9.3 X 10-'^ 
7.5 X 10-5 

1.3 X 10-5 
548 

7.3 X 10-5 

1.7 X 10-5 

1.8 X 10-5 

6.9 X 10-'^ 
138 

1.5 X 10-5 

2.2 X 10-5 
3.1 X 10-5 

4.3 X 10-5 
3900 

2.0 X 10-5 
2.8 X 10-5 
4.5 X 10-5 
7.7 X 10-5 
3400 

2.7 X 10-5 
2.5 X 10-5 
4.2 X 10-5 
1.1 X 10-5 
1900 


Table 4: Discovery prospects for the three benchmark points. £ denotes the integrated 
Inminosity reqnired for a detection in the given channel. 


these tend to increase with increasing Hi mass, we note that for H 2 —>■ in point 2 and 

3 the VBF channel is somewhat better. 

Comparing the detection prospects in tablewith the reach shown in the corresponding 
plots, it is clear that the additional constraints on the kinematics dne to the assnmption 
abont the mass spectrnm are important in improving sensitivity, thongh in some cases the 
differences are negligible. The invariant mass cnts are especially important for WH as the 
ttbb and tt backgronnds increase signihcantly with increasing overall invariant mass, hence 
rendering a (relatively low) cnt on over all invariant mass very effective. 

To illnstrate more clearly the signal and backgronnds in the WH channels we plot these 
as fnnctions of the invariant mass of the 6-jet pair^ in hgnre [TT] for Hi —)• AiAi in 
point 1 (hgnre [T^a)) and for H 2 —?■ in point 3 (hgnre [TT|(b) ). Note that, dne to the 

cnts in the analysis, the distribntions are restricted to ±15 GeV and that both plots 
nse single 6-jet analysis only. In hgnre [T^a) it is clear that a somewhat narrower cnt wonld 
signihcantly rednce the backgronnd from ttbb withont affecting the signal signihcantly. 
However, it is important to note that there are signihcant statistical nncertainties in the 
backgronnds (partially hidden by the smoothening) and this, together with the fact that 
the signal already dominates and that we do not inclnde detector resolntion, means that 

^Included here are the two invariant masses of the combination of fe-jet pairs that have the smallest 
difference in invariants mass of all the combinations of 6-jet pairs where both pairs have invariant mass 
within mAi ± 15 GeV. 
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Figure 11: Signal and backgrounds for two of our benchmark points in the WH channel. 
The backgrounds have been smoothened for visual clarity. 


it is not clear how beneficial such a cut would really be. From figure [T^b) we see that 
for heavier masses the invariant mass peak is smeared out towards smaller masses, this is 
mostly an effect of the difficulty in finding the “right” 6-jet pairs; for lighter ^is the “wrong” 
pairing gives too high invariant masses to be included. 

9 Prospects at Higgs factories 

Since we have seen that the conditions for detecting light pseudoscalars at the LHC via 
VBF and HS are rather challenging and it may require several years to do so, it is worth 
investigating whether a Higgs factory would do better in this respect. To this end, we 
estimate the sensitivity of a 240 GeV e'’'e“ collider that produces CP-even Higgses primarily 
through the ZH channel. Since the hadronic background here is much smaller than at the 
LHC, we do not require any 6-tagging but just look for a leptonically decaying Z-boson 
(into electrons/muons) accompanied by exactly 4 jets (each with px > 15 GeV). 

In all other aspects our analysis is identical to what was done for the LHC, except that 
we do not include jet substructure techniques at this point. In fact, we only study chains 
starting with iLsM for which we assume a production cross section of (T(e^e“ —>■ H^y^Z) = 
200 fb. The resulting sensitivity is shown in figure It is clear that a Higgs factory would 
be much more suitable for our channels, as sensitivity to the NMSSM dynamics will be 
established already in the first year of operation. Furthermore, the very low background 
also means that profiling the events in terms of underlying mass spectra and H^mWW 
coupling strength would be an easier task than at the LHC. 
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Figure 12: Reach of a 240 GeV e'’'e“ collider in H 2 —)• for H 2 = F^sm in the ZH 

channel. All sensitivity curves assumes mH 2 — 125 GeV and looks for four jets with total 
invariant mass of 125 GeV in addition to the leptonically decaying Z. The colour code for 
the points is the same as in 



10 Conclusions 

In contrast to the more constrained MSSM, the NMSSM allows for the existence of very 
light Higgs scalars as well as pseudoscalars. Therefore, the discovery of, in particular, a 
light pseudoscalar Higgs state would not just prove the existence of physics beyond the SM 
but would also be inconsistent with minimal supersymmetry. 

Since the NMSSM also accommodates the 125 GeV Higgs more naturally than the 
MSSM, it is well worth a study. Although the most promising LHG discovery channels 
of the aforementioned light pseudoscalar state are most likely based on GF production 
of heavier scalars that subsequently decay to AiAi or AiZ, we demonstrated here that 
also VBF and HS production of the heavier scalars can be accessible. Hence, these two 
additional production modes can be exploited to study couplings not accessible in GF, such 
as those of the heavy scalars to SM gauge bosons. Especially interesting are the channels 
starting with Hs'- as these can have BR(Ffj —?■ AiA\) close to 1, the channels of this paper 
might be our only chance of measuring the couplings of these scalars to gauge bosons. 

In these channels, the signal rates are substantially lower than in the case of GF, but 
the same is true for the backgrounds. Due to the nature of the couplings involved, the 
only decay chains of interest here are Hi ^2 AiAi. For VBF production of Hi ^2 the most 
promising of our final states is 262r (in addition to the two forward/backward jets), which 
may allow detection, although only at 3000/fb. 

For HS production, the background is even smaller, especially for the ZH mode. How¬ 
ever, this channel also has a very small cross section and hence the signal will be very hard 
to extract. As the event rates are indeed significantly smaller compared to VBF and GF, 
here it is most beneficial to employ the final state with the highest rate, i.e., 46. 

Although still featuring a relatively low signal rate, WH production shows significantly 
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better prospects than ZH. Despite significant backgrounds from ttbb and tt, the higher 
signal rates make this the most promising channel studied in this paper, at least for relatively 
light initial scalars. However, given the invariant mass structure of the main backgrounds, 
the signal tends to be overwhelmed unless one can cut on the four-6 invariant mass and 
the enforced mass window needs to be relatively narrow. Therefore, the prospects for this 
channel diminish as the mass of the initially produced scalar increases, rendering VBF the 
most promising channel for heavier scalars. 

In addition to general scans for sensitivity reach in parameter space, we have performed 
more detailed studies of three representative benchmark points. Here, we assumed knowl¬ 
edge of the masses of the produced scalar as well as the pseudoscalar (e.g., as measured in 
GF production) and used this information to constrain the kinematics. Especially for HS, 
this can dramatically improve the sensitivity. 

Finally, we have briefly tested how the LHC prospects compare to those of a future 
e'*'e“ machine running at the threshold of ZH production, the golden (Bjorken in this 
case) channel of such a Higgs factory. Needless to say, sensitivity to decays 

is rather prompt herein. Further, the clean environment of an e^e“ machine in terms of 
very limited backgrounds will enable precise diagnostic of the NMSSM dynamics involved. 
It remain to be seen whether commissioning of such a leptonic machine will occur on a 
timescale that will make it competitive with a potential LHC stage at high (instantaneous) 
luminosity, the so-called Super-LHC (SLHC) option [92], currently being discussed as a 
probable upgrade of the CERN collider past Run H. 
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